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We investigate the flow properties of a two-dimensional aqueous foam submitted to a quasistatic 
shear in a Couette geometry. A strong localization of the flow (shear banding) at the edge of the 
moving wall is evidenced, characterized by an exponential decay of the average tangential velocity. 
Moreover, the analysis of the rapid velocity fluctuations reveals self-similar dynamical structures 
consisting of clusters of bubbles rolling as rigid bodies. To relate the instantaneous (elastic) and 
time-averaged (plastic) components of the strain, we develop a stochastic model where irreversible 
rearrangements are activated by local stress fluctuations originating from the rubbing of the wall. 
This model gives a complete description of our observations and is also consistent with data obtained 
on granular shear bands by other groups. 
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Cellular materials, such as foams, concentrated emul- 
sions, slurries, or granular materials, exhibit Theological 
properties that cannot be understood within the scope 
of standard solid or liquid mechanics [EL 0, pi . For such 
systems, thermal energies are orders of magnitude lower 
than the typical energy required to relax the structural 
arrangements of their components; under small forces, 
the material remains trapped in a metastable configura- 
tion and exhibits a solid-like behavior. When submitted 
to a large enough stress however, it can be driven through 
a sequence of new metastable configurations, giving rise 
to a macroscopic flow. But the resulting flow field may 
still differ a lot from what would be expected for a molec- 
ular liquid. 

Dry sand slowly flowing down an hourglass provides 
a simple example of such abnormal flow behaviors: the 
flow splits into a plug-like central region and a strongly 
sheared thin layer at the wall - a few particles wide - 
where most of the dissipative process occurs 0. This 
spontaneous localisation of the strain in narrow regions 
of the material (the so-called shear bands) can be ob- 
served in many other situations such as shear, surface, 
or convective flows for instance §|, @, || ■ Shear band- 
ing actually controls most of the practical situations one 
has to face in soil mechanics and industrial handling of 
grains, and is also relevant to pyroclastic flows in geology 
(for a review on granular matter see ||). This question 
has recently received a lot of attention from physicists, 
both theoretically and experimentally [jl0| [llj , but a clear 
picture has not emerged yet. 

By contrast, the possibility of shear banding in foams 
has been mostly ignored in the literature, and numerical 
or theoretical studies usually assume shear flows in foams 
to be uniform |2, 12 . The assumption that shear banding 
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is unique to granular matter can be misleading because 
it suggests that some peculiar aspects of granular flows, 
such as solid friction, particle rotation or dilatancy, are 
required to derive a shear band model. 

In this Letter, we report the formation of shear bands 
in aqueous foams. We believe that foams may shed light 
on the dynamics of granular systems by evidencing the 
minimal set of ingredients needed to get shear banding. 
To that extent, foams constitute a much simpler model 
than granular systems since the basic bubble/bubble in- 
teractions which control the mechanical properties of the 
material are well known: elastic (stored) energy is related 
to an increase of the total interfacial area when the bub- 
bles are distorted whereas dissipated energy is associated 
with neighbors swapping events (Ti processes) inducing 
flows in the liquid films and vertices (for a review on 
foams, see Q). 

In order to probe the microdynamics of the foam, one 
needs to track the trajectory of each bubble during shear- 
ing. Since 3-D foams are inherently diffusive to light, we 
used a 2-D model foam - a monolayer of bubbles - sub- 
mitted to a continuous slow shear in a Couette geome- 
try. The setup was composed of an inner shearing wheel 
and an outer ring (of respective radius Rq = 71 mm and 
i?i = 122 mm) confined between two transparent plates 
separated by a 2 mm gap. To produce the foam, the 
cell was first hold vertically and partially filled with a 
controlled volume of soap/water solution. Bubbles were 
formed by blowing nitrogen gas through two small injec- 
tion holes at different flow rates until the resulting foam 
reached the top of the cell. Once set horizontally, the 
foam rapidly attained a uniform wetness characterized 
by its liquid fraction 0.01 < < 0.3 (Fig. 1). This foam- 
ing procedure was chosen because it produces bidisperse 
disordered foams and therefore eliminates crystallisation. 
The mean diameters within each of the two populations 
of bubbles were of the order of 2 and 2.7 mm, with a mean 
deviation of 0.2 mm. These bubbles were large enough 
compared to the gap height so that they would not over- 
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FIG. 1: Close-up frames of dry and wet 2-D foams under con- 
tinuous shear. The shearing wheel appears in black. Notice 
the teeth shaped edge which eliminates slippage of the first 
row of bubbles. Arrows indicate the shearing direction. (A) 
A dry bidisperse foam (<j> = 0.05), showing deformed polygo- 
nal cells. (B) A wet bidisperse foam (<f> = 0.20). The bubbles 
are circular and undeformed. (C) System of coordinates used 
to analyze the flow field. 
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FIG. 2: Average tangential velocity profiles for different liquid 
fraction (♦:<?(> = 0.05, A: <j> = 0.12, •: <j> = 0.20), showing 
exponential decays with a width A (</!>). The small plateau 
for < r < 1 corresponds to the first row of bubbles being 
attached to the shearing wheel. Inset : dependence of A on 
liquid fractions cj> (the line is just a guide for the eyes). Two 
plateaus can be distinguished on either side of <j>* ~ 0.12, 
which separate deformed and undeformed bubbles regimes. 



lap and form a truly two-dimensional foam. To define 
a bubble scale, we measured the mean distance d be- 
tween first neighbors in the foam. In all experiments, 
d lay between 2.1 and 2.5 mm so that the gap between 
the wheel and the ring could accommodate from 20 to 
25 rows of bubbles. The distance d was evaluated several 
times during the experiment and found to be almost con- 
stant. Coarsening would eventually lead to a growth of 
the biggest bubbles at the expense of the smallest ones, 
but over a longer time. We also checked the absence of 
shear induced size segregation that might have occurred 
during the experiment. 

Shearing was induced by rotating the inner wheel at 
constant velocity V w heei using a stepper motor. To avoid 
slippage at the wheel and the ring, their sides were teeth 
shaped so that the first and last rows of bubbles would 
remain irreversibly attached to the walls. To eliminate 
transient effects, we ran the experiment a full round be- 
fore taking data. The motion of 1000 to 1500 bubbles 
was then recorded using a CCD digital camera positioned 
over the setup. In a typical experiment, 3000 images were 
taken corresponding to a total displacement of 600 d of 
the wheel edge. The apparent centers of mass of the bub- 
bles were subsequently tracked by image analysis (IDL 
software). To reduce the effect of the viscous friction be- 
tween the bubbles and the confining plates, we restricted 
our study to quasistatic flows. We focused on average ve- 
locity measurements as a probe of shear rate dependence: 
we found that in the range < V w heei < 0.7 mm. s -1 , the 
velocity profiles were similar apart from an overall scale 
factor. All experiments were performed in the quasistatic 
regime at V w heei = 0.25 mm. s -1 . In the following, all 
velocities and distances are normalized by V w heei and d 



respectively. We note ujq = V w heei /d the characteristic 
frequency of the shear. 

Figure 2(A) shows the decay of the average tangen- 
tial velocity (vg(r)) with the distance r to the shearing 
wheel for different liquid fractions <f> (see Figure 1(C) for 
variables definition). Averaging was performed over the 
tangential coordinate rg and time i, yielding smooth and 
reproducible profiles, although the instantaneous flow is 
strongly intermittent. The reduced velocity is found to 
approach 1 at r — » confirming the absence of slip at 
the edge of the wheel. At larger r, the profiles exhibit an 
exponential decay : 

(v e (r)) ~ exp(-r/A) (1) 

with a width A depending on <f>. The curve A versus 
<fi, presented on Figure 2(B), shows two plateaus at low 
and high volume fraction. The transition between these 
two regimes occurs around <jf = 0.12 which qualitatively 
marks the limit between dry foams with polygonal bub- 
bles for tf> < <f>* and wet foams with undeformed bub- 
bles for (/)>(/)*. In both cases, the rapid decay of the 
mean velocity over a few bubble diameters, establishes 
the existence of shear banding in foams. The exponential 
shape of the velocity profile, observed in all experiments, 
appears as a robust feature which was also observed in 
comparable experiments performed on 2-D granular ma- 
terials § § @. 

Beyond these time averaged profiles, the present setup 
allows measurements of the short timescale fluctuations 
of the bubbles velocities. A mere observation of the video 
sequences reveals brief oscillations of clusters of bubbles 
of various radial extension, rotating together as rigid bod- 
ies as shown in Figure 3. These dynamical structures 
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FIG. 3: Video frame of the foam with the position of the 
bubbles centers and their trajectories over the last 20 seconds. 
This time period corresponds to a total displacement of one 
bubble diameter for the first row (or equivalently for the inner 
disc edge). The dots size and lines color reflect the total 
distance travelled by the bubbles revealing a large rotating 
cluster. 



are ephemeral and disappear after the wheel edge has 
moved by roughly one bubble diameter (this was checked 
by measuring time correlations of the velocity which we 
found decay to in a time of the order of 1/ujq). To 
quantitatively probe these coherent moves, we studied 
the spatial correlations of the instantaneous velocity field. 
We focused on the radial component v r which has a zero 
time average and therefore gives a better signal to noise 
ratio (qualitatively, similar results were found when using 
vg — (vg) instead of v r ). Figure 4(A) shows the correlation 
function g r ( Arg) = (v r (r, rg).v r (r, rg + Arg)) / (vl(r)) for 
different values of r from 1 to 10, at a volume fraction 
4> = 0.20. Regardless of r, g r decreases with Arg from 1 
to a negative value then slowly relaxes to 0. The length 
£(r) for which g r reaches defines a typical correlation 
length of the velocity field at a distance r. In Figure 4(B), 
rg has been rescaled by r. All the curves then collapse on 
a single one, which demonstrates a linear increase of £(r) 
with the radial distance r: £(r) = ar. Motivated by the 
observation of oscillating clusters, we modeled the flow 
field as a superimposition of rotating blocks of bubbles of 
various sizes (see figure caption for details) which allowed 
us to obtain a good fit of the master curve (Figure 4(B)). 
Similar results were obtained at all volume fractions, but 
the coefficient a was found to decrease with 4>. 

It is crucial to note that the velocity field which yields 
the measured correlation functions mainly corresponds 
to reversible movements of the bubbles centers, and thus 
probes the elastic deformation of the foa m rather than 
the plastic flow. The quantity \J < v% >/ < vg > pro- 
vides a good estimate of the ratio of reversible to irre- 
versible moves of the bubbles occurring upon shearing. 
This quantity is larger than 1 beyond the first attached 
row of bubbles, and gets larger than 10 beyond the fifth 
row. The correlation measurements thus reveal that the 
instantaneous stress field is spatially correlated. 

This peculiar characteristic of the foam deformation 
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FIG. 4: (A): Spatial correlations of the radial velocity for 
different radial distances r from 1 to 10 (<f> = 0.15). (B): All 
correlations can be collapsed on a single curve when plotted 
versus the tangential distance rescaled by r. To fit this master 
curve, we assume that the radial velocity v r (r, re) at different 
r comes from the rotation of clusters of mean lateral extension 
£(r) proportional to r (£(r) = ar). The radial velocity field 
within a cluster of size /i is idealized by a two step function 
(v r = evo for — fi < re < and v r — — evo for < rg < /.t with 
e = ±1). We postulate a statistical distribution of /j, around 
£(r) by assuming that the fraction of clusters of size larger 



than p is exp 



The resulting correlation function 
ar I c ^ I ar I ■ 



(dotted line) then reads: g r (Ar e ) = ( 1 - ^) exp 



field can actually be understood under the scope of lin- 
ear elasticity, by simply modeling the foam as an isotropic 
elastic medium. During the initial loading, a uniform 
mean stress a develops in the material (we neglect the ra- 
dial geometry since the wheel radius is much larger than 
the shear band width). In the steady state, this uniform 
stress persists in average but is locally modulated by a 
fluctuating stress field Aa{t) of mean value associated 
with the continuous rubbing of the wheel teeth. At each 
location (r = 0, rg — i) on the wheel edge, the foam is in- 
deed submitted to a localized perturbative stress Act.; (t) 
of variance sq = (Actq) varying at a frequency uiq which 
elastically propagates into the material. These multiple 
noise sources add up to produce a stress fluctuation at a 
position (r, rg ) of amplitude (see for instance fia ) : 



Aa^t) 



\Jr 2 + (rg - i) 2 



(2) 



Assuming the noise sources to be uncorrelated 
({A<Ji(t)A<Tj(t)) = for i ^ j), the resulting stress co- 
herence length, at a distance r, takes the form £(r) = ar, 
in agreement with our experimental findings. We inter- 
pret the different observed values for a as a signature of 
the anisotropy of the foam due to its initial loading. This 
is consistent with the observation of large values of a for 
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the driest foams where the largest uniform deformation is 
first produced. From Eq. ^, we are also able to compute 
the variance of the fluctuating stress Ao> (t) as: 



< Aa 2 r >= s(r) ~ ^ 
ar 



(3) 



At this point, we wish to relate the fluctuating stress 
field Aoy(i) to the measured average flow profile pre- 
sented before, by taking into account the plastic property 
of the foam. Our approach mostly follows a model pro- 
posed by Pouliquen and Gutfraind H based on Eyring's 
activated process theory 16 to describe chute flows of 
granular materials. In the present description, the vari- 
ance of the local stress fluctuation s(r) plays the role of a 
temperature allowing plastic flow to occur. The moving 
boundary acts as a "hot wall", exciting internal defor- 
mation modes of the material in the form of self-similar 
rotating clusters. When the fluctuating stress overcomes 
a certain yield stress a y , the structure plastically yields. 
The yielding rate in the material by unit of time and 
space thus writes : u) = loqP{g > a y ) where P(<j) is the 
density probability of stress. The stress at a distance r 
is a sum of ~ r random variables (see Eq. ||) so that 
P(<r(r)) is a gaussian distribution centered on a of vari- 
ance sir). Using Eq. [| the yielding rate at a distance r 
reads: 



lo{t) = cj Q P(a(r) > a y ) = lu 



I -<TfJ- 



with 
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(4) 



(5) 



Each failure increments the average velocity gradient by 
1 in reduced unit, so that the constitutive equation for 

the flow writes — ~ ~ — u(r)/wo with the boundary 
conditions (ve(0)} = 1 and (vg(oo)) = 0. A very good 



approximate function to the integral of w[r) is given by 
a pure exponential so that (vg(r)) ~ exp(— r/A). 

This model gives the right form of the tangential ve- 
locity decay. An outcome of the model is contained in 
Eq. H which relates the a coefficient to the plastic flow 
decay length A in the foam. Experimentally, we found 
aX = 1.02 ± 0.08 (A ranging from 1 to 2), which would 
suggest that so/(a y — a) 2 is independent of the liquid 
fraction eb. This quantity is indeed self-adjusted via a in 
the transient regime to allow enough rearrangements to 
occur. 

In conclusion, we have observed shear banding in dry 
and wet 2-D foams under continuous slow shear and we 
have probed the associated elastic deformations of the 
foam, characterized by brief, collective oscillations of self- 
similar blocks of bubbles. We have developed a stochastic 
model which relates the plastic flow to the stress fluctua- 
tions experienced by the foam. The main characteristics 
of the flow (rapid decay of the average velocity over a 
few bubbles, large velocity fluctuations) are very similar 
to what is commonly observed in granular systems, sug- 
gesting that the proposed mechanism could remain valid 
for granular systems. As already mentioned, the pre- 
dicted exponential velocity decay has been observed in 
various 2-D granular shear bands j|, |[ FS. Moreover, the 
velocity profile in 3-D has been shown to obey a gaussian 
decay in the limit of disordered and non spherical grains 
Q . This functional form for the velocity profile immedi- 
ately follows from the modification of Eq. || in 3-D which 
then writes s(r) ~ ^° 2 yielding a gaussian decay for 
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